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Abstract 

Q^' We propose that dark matter is dominantly comprised of atomic bomid states. We build a simple 

' \ model and map the parameter space that results in the early universe formation of hydrogen-like dark 

atoms. We find that atomic dark matter has interesting implications for cosmology as well as direct 
\ detection: Protohalo formation can be suppressed below Mproto ~ 10'^ — lO^M© for weak scale dark 

matter due to Ion-Radiation interactions in the dark sector. Moreover, weak-scale dark atoms can ac- 
commodate hyperfine splittings of order 100 keV, consistent with the inelastic dark matter interpretation 
^ ■ of the DAMA data while naturally evading direct detection bounds. 

m 
in 

■ 1 Introduction 

O 

I Cosmological observations suggest that dark matter comprises more than 80% of the matter in the universe 
[H [2j . Much of the effort to explain the origin of dark matter has focused on minimal solutions in which 
, dark matter consists of a single particle species, the most popular being the neutralino in variants of the 
supersymmetric standard model. Such dark matter models include the compelling feature that weak-scale 
physics ~ weak-scale mass and weak-force coupling strength - can naturally generate dark matter with the 
correct cosmological abundance. Dark matter in this broad class is described as weakly interacting massive 



X 

H 

. ; particles (WIMPs). 



However, conflicts do exist between WIMP models and observational data. Simulations of WIMP dark 
matter predict significantly more small-scale structure than current observations suggest [3l[l]. In addition, 
the direct detection experiment, DAMA [5], sees a positive signal with great significance (8ct), yet when 
interpreted as a standard WIMP, other experiments such as CDMS [6] and XENONIO [7], completely rule 
out the same parameter space. Finally, measured cosmic ray spectra may suggest a new primary source 
for electrons and positrons in our galaxy and potentially evidence for dark-matter annihilation; however, 
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the standard neutralino candidate is unable to fit this data [HI [9l [TOl [TT]. 

These issues suggests compelling reasons to explore dark matter models beyond the minimal candidate. 
In addition, the dark matter sector (or 'dark sector') may be rich with complexity and may feature unan- 
ticipated dynamics. In fact, the dark matter may even interact via a long-range force - a massless gauge 
boson - which is still allowed by the bounds on the number of relativistic degrees of freedom during big 
bang nucleosynthesis [12] . 

In this paper we propose a dark sector charged under a hidden U(l) gauge symmetry. We assume two 
species of fermions, a 'dark proton' and a 'dark electron', and that the dark matter abundance comes from 
a matter-anti-matter asymmetryj^ We shall see that in interesting parts of parameter space, the bulk of 
the dark matter exists in atomic bound states. The Lagrangian is 

^darfc ='^p(^ + "lp)^'p+'^e(-P + "2e)^'e (1) 

where P = igl + gQ^ and Q = ±1 for ^'p and ^'e respectively. In what follows we use the convention 
nip > me without loss of generality. We show (Section 2) that for parts of parameter space, recombina- 
tion in the dark sector occurs efficiently, and we discuss the bounds from and implications for structure 
formation. We then add interactions which allow for direct detection in a way that mimics inelastic dark 
matter [15] and show that there exist parts of parameter space which can explain the DAMA signal, while 
avoiding constraints from other direct detection experiments (Section 3). Finally, in Section 4 we discuss, 
in a cursory way, other phenomena potentially related to atomic dark matter. 

A number of ideas related to this work have appeared in the literature. For example, the idea of U(l) 
charged dark matter has appeared in [W\ \17\. the idea of composite dark matter in |18j . and that of mirror 
dark matter in [191 120j. To our knowledge, this is the first work to explore the generic parameter space for 
viable atomic dark matter. 



2 Cosmology 

Introducing a new hidden U(l) has interesting cosmological implications. Our interests lie in the parameter 
space that affords atomic systems. The existence of standard model (SM) atomic hydrogen states in the 
early universe requires an asymmetry between particles and antiparticles; dark atoms are no different. 
We assume that there is a 'dark asymmetry' akin to the baryon asymmetry in the SM, and that the 
dark asymmetry is such that the universe is net charge neutral, rie = repH The existence of dark atoms 
implies that dark matter is coupled to dark radiation until the universe cools beyond the binding energy 
of hydrogen 

B = \alm, (2) 

where ud is the dark fine structure constant and /xh = {memp)/{me + mp) is the reduced mass of dark 
hydrogen. This has potentially interesting implications for structure formation because interactions in the 

^ Some models that use the matter-anti-matter asymmetry to generate the correct dark matter abundance exist [131 114] , 
but we do not explore them here. 

^Unless otherwise noted, e, p, and H refer to the dark electron, dark proton, and dark hydrogen, respectively. 
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dark sector can decouple much later than in a conventional CDM WIMP model. Observations of satellite 
galaxies seem to favor some mechanism to damp the growth of small scale structure in dark matter [211 [22] , 
which, as discussed below, can be provided by atomic dark matter. 



2.1 Dark Recombination and Halo Constraints 



One of the most interesting features of the model is the presence of both neutral and ionized dark mat- 
ter components. The fractional ionization, Xg, plays an important role in the cosmic evolution of the 
dark matter. At early times, X^. affects the decoupling temperature of dark matter and dark radiation, 
which impacts small-scale structure formation of dark matter. At late times, bounds on dark matter self- 
interactions constrain X^. because the dark matter ions interact through a long range force. The residual 
ionization fraction in the dark sector is governed by neutral atom formation in analogy with SM hydrogen 
recombination [23]. In the following, we follow the notation of Ref. 



The residual ionization fraction is found by solving the Boltzmann equation for the free dark electron 
fraction, 



X^ 



Up 



He + nu 



(3) 



The evolution of X^ depends on the Hubble rate, H, and the rate for e + p ^ H + 7. We can write the 
thermally-averaged recombination cross section using the dimensionless variable X rJ-^ ctS 



(av) = ^ 



/277r /ifi 

where ^ = 0.448 is a best-fit numerical coefficient \25\ 126). The equation governing X^ can be written as 



(4) 



dX. 



dx Hx 



(1 - X^fp - XlnoMiov) 



(5) 



where 



/3 = {<yv) 



Brue 

2-KX 



3/2 



(6) 



As discussed in [231 125] . recombination into the n = 2 state completely dominates the evolution of X^. 
This is accounted for through the factor C in Eq. ([5|) which represents the fraction of n = 2 states that 
produce a net gain in the number of ground state hydrogen atoms. This is not unity because the thermal 
bath can ionize the n = 2 state before it decays. Thus, C is the ratio of the (n = 2 — > n = 1) decay rate 
to the sum of this decay rate plus the ionization rate (see [23] for a detailed discussion) 



C 



K + A27 

A„ + + /3(2) ' 



(7) 
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Figure 1: The allowed parameter space in mg for a given nip and as a function of the residual ionization 
fraction, Xe- Atom dark matter is viable in the colored regions, which correspond to 10^^ < < 10^^ 
(red circles), 10^'^ < Xe < 10^^ (blue triangles), 10"'^ < X^ < W^^ (magenta boxes) and Xe < 10^"^ 
(green stars). The striped region is ruled out by Eq. [21] with k = 3 and this region extends to the dashed 
black line for n = 10. The black-crossed region is ruled out because Xe > 10%. 



The rates are given by 

A - (9) 

- (8vr)2(l -^e)n^M ^ ' 

A2, = 3.8^^ (10) 

where the rate for a Lyman-a photon to redshift such that it cannot excite n = 1 ^ n = 2 and A2'y 
is the two photon decay rate of 2s Is, which has been taken from Ref. [27J. We find that X^ varies 
from 1 — 10"^'^ throughout the parameter space od S [10~^,0.3] , rrie € [0.01 GeV, mp], nip G [me,3TeV]. 
Self-interactions, as discussed below, rule out some of this parameter space. A few representive planes are 
plotted in Figure [TJ 

Bounds on the atomic parameter space and Xe can be derived from observations of the Bullet Cluster and 
halo profiles. The bounds are derived through the momentum transfer cross sectioiJl 

al, = I dl^(l-cos0)^ (11) 



■^Note that amt reduces to the total elastic cross section for hard-sphere s-wave scattering which is typical of WIMP dark 
matter models. 



where the index i runs over the three types of self-interactions present in our dark sector: Hydrogen- 
Hydrogen, Ion-Hydrogen, and Ion-Ion. The last process is described by Coulomb scattering, but since we 
want to study a dominantly atomic dark sector, the Ion-Ion cross section is the least relevant to our model 
and we do not discuss it further. Naively, the first two cross sections are bounded by geometric values 



where oq = l/(a /^h) is the Bohr radius of dark hydrogen. However, this naive guess is inadequate. At low 
energies {k <^ 1) both of these processes can be described by scattering from a central potential 



V oc r"" (13) 

with n = 6 and 4, respectively. In this case one finds that the cross sections are velocity-independent 
constants and enhanced over the geometric estimate [28]. We are generally interested in a wide range 
of kao, thus these results are not strictly applicable however these cross sections are slowly decreasing 
functions of the relative velocity ^28j. A conservative estimate of the cross sections is given by 



I < (-0)^ (14) 

with 3 < K < 10. The values of k have been inferred from general quantum mechanical scattering [29\ [5U] 
and detailed computations of SM hydrogen scattering [311 [32l [33] . 

With the relevant cross sections in hand, we can use observations of the Bullet Cluster [Ml [35] as a guide 
for the present day maximum value of Xe- Measurements of the mass-to-light ratio and the radius of 
the sub-cluster suggest that the sub-cluster could have lost no more than Fobs = 20 — 30% of its initial 
mass. Following the analysis in [M], the number of scattering centers that a single dark matter particle 
encounters as it passes through the target cluster in the case of one species is 

r = -Es; (15) 
m 

this quantity is often referred to as the scattering depth. The parameter is the surface mass density of 
the sub-cluster defined as 

'^piz)dz, (16) 



where p{z) is the sub-cluster's volume mass density and R is the radius of the sub-cluster. For multiple 
species, with species i in the sub-cluster having a mass density Hj, and mass rrii, scattering off of species j 
in the target cluster, we have 



m 



aijfj (no sum) (17) 
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where aij is the cross section for i and j to interact and fj is the number fraction of species j in the target 
cluster. Equation [T7] can be rewritten in terms the of total observed surface mass density, Hj^ = ^ Hj, as 

i 

Tij = ^T^^=Jz O'ij/j- (18) 

2^ JfcWlfc 
k 

We make the simplifying and conservative assumption that all of the ions in the sub-cluster are scattered 
out of the sub-cluster. In this case, the mass fraction lost through H — j scattering is bounded by 

Fn = Fobs - /p^ - fe— = Fots - Xe. (19) 
mn niH 

The mass fraction actually lost from the sub-cluster, using cross sections as parameterized in dHD and (fTSj) . 
is 

^ = T.^Hj = ^{l- X^)A7T^al (20) 
Demanding that A < Fn we have the bound 

a ^ (21) 



mn mn St(1-X|)' 

Plugging in the values = 0.2-0.3cm^^g and Fobs = 0.2 gives a constraint on the atomic parameter 
space [M] 

OlWlGeVy/lOOOeVX 

an/ V / V run ) ^-^i 

Thus, we find that and cr/m-H are bounded simultaneously. From our conservative (and representa- 
tive) assumption about the Ion-Ion cross section, is bounded to be less than 10% — 20% regardless of 
a/mn- For very small X^., the usual CDM WIMP bounds on cr/mn are applicable; our estimate yields 
cr/nin < Icm^/g, which is slightly larger than the detailed simulations of |35j. Figure [1] shows some of 
the allowed parameter space in the an — rne plane for a few atomic masses ranging between lOGeV and 
ITeV. 



Previous considerations of a hidden U{1) |16tll7j have concluded that soft scattering of charged dark matter 
can drastically affect halo formation and thereby rule out large swaths of parameter space. This result 
follows from the soft singularity in the Rutherford scattering rate which, when integrated over galactic time 
scales, can lead to significant energy transfer between charged particles. This effect tends to smooth out the 
core of the dark matter distribution. Application of these results excludes all of the parameter space shown 
in Figuredl However, since Hydrogen- Hydrogen scattering is well modelled by hard-sphere scattering in 
the majority of the considered parameter space, these bounds are not applicable to atomic dark matter. 
The relevant bounds from halo formation considerations are O.lcm^/g < cr/m^ < Icm^/g, which do not 
signficantly change our conclusions [36]. Atomic dark matter provides a dynamical mechanism to shut off 
the naive long range effects of a hidden U(l). 
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2.2 Protohalos and Radiation Damping 

Certain values of the residual ionization fraction may enable the ionized dark matter to smooth out structure 
on small scales. The parameter space of atomic dark matter allowed by constrains derived in Section 12.11 
predict that kinetic decoupling of the dark matter fluid and dark radiation occurs during the radiation 
dominated epoch. The dark radiation damps the growth of structure until it decouples from the dark 
matter [37] . The decoupling temperature is given by 

T{Tdec) = H{Tdec) (23) 

where 

r(r) = n^ar = X^udmctt (24) 

H{T) = ^ 5.5— (25) 



and (Tt is the Thomson cross section, cjt = 8/37r(Q£)/me)^. Using ([23l) - (p5]) the decoupling temperature 
is given bjfl 

A full analysis of the power spectrum is left to future work (see |16 1 I37 1 [38] for discussions of similar effects), 
however, a rough estimate of the comoving wavenumber that is damped due to the radiation interaction is 

kdamp ^ (27) 
Vdec 

where rj^ec is the conformal time of decoupling. This can be written in terms of today's temperature, Tq, 
and the decoupling temperature Tdec 

r\dec = 2W(1 + Zdec) ^ ^ ^ ■ (28) 

The mass scale which survives damping, and which eventually characterizes the first dark matter clumps, 
can be written as 

Mgrow > — ( -;- ) ^DMPcrit 

where pcrit is today's value and the largest value of Mgrow corresponds to kdec = Vdec- "^^^ large range 
of Mgrow occurs because we have allowed for an order of magnitude error in kdec- Thus, we find that 
atomic dark matter can have much less power on small scales compared to a conventional CDM WIMP if 
there is a sizeable fraction of free electrons. In particular, weak-scale dark atoms can be consistent with the 
observed populations of intermediate mass dark halos in the Milky Way [1] . We emphasize that simulations 
are necessary to evaluate the detailed power spectrum and satellite popluations predicted by atomic dark 
matter. However, it is clear that even the simplest atomic dark matter system can significantly impact 
structure formation. 




Most of the allowed parameter space requires B < Aqcd, thus we take g, = 12. 
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3 Direct Detection 



Atomic dark matter, as thus far considered, is secluded from the standard model. While the cosmology of 
atomic dark matter is interesting in its own right, it naturally lends itself to inelastic scattering because 
of energy level quantization. This offers an exciting possible explanation of the DAMA data [5]. 

The unperturbed energy levels of hydrogen are 

^» - (30) 

One might hope that the DAMA scale - 0(100 keV) - could be generated by energy differences between 
levels with different principle quantum numbers. Generically the rate of elastic scattering will be greater 
than that of inelastic scattering. However, predominantly inelastic scattering could be enforced by setting 
irip = rrie = In this case, the first Born term for elastic scattering vanishes. Unfortunately, efficient 

recombination in such a scenario forces one to consider mn ~ GeV, which is too small to account for the 
recoil energies measured by DAMA. Nevertheless, atoms have a rich structure and the allowed parameter 
space for viable recombination naturally leads to hyperfine splittings on the order of 100 keV for weak-scale 
hydrogen. The hyperfine splitting is given by 

Ehf = i^gegpOiDf^H — (31) 

where ge,gp are the gyromagnetic ratios of the dark electron and dark proton, which we take to be equal 
to two. 

Exploiting this scale requires a scattering process which is dependent on the spins of the dark atom's 
constituents. This can be accomplished with a broken U{l)x which is axially coupled to the dark matter 
and mixed with the standard model hypercharge as in jl8j 

^mix — i Bp^p. (32) 

Having an axial coupling in the dark sector and a vector coupling to the standard model will ensure that 
the dominant scattering process changes the dark atom spin state by one unit. After integrating out the 
Z boson and diagonalizing the gauge kinetic terms, the Lagrangian is 

= ^^SIVI + CdM + ^Dark Gauge 

^ g 

CdM = ^'p(i^-g5 75^+5/ + "T'p)^p + ^e(«^ + 5'5 75^-5'4 + "ie)^e Y fiJr) 

J-DarkGauge — ^ i^-i^ 4 ^"^ \^^w^EM~^^^w\^^^j 'Jzj^fJ-' 2 ' ^ 

The parameters and are the cosine and sine of the weak mixing angle. The dark current Jj^ is 

Jd = -9n "^pT^Ts^p + 55 '^eT'^Ts^e, (34) 
and J^j^j and are the standard model electromagnetic and weak neutral currents, respectively. 



8 



The calculation of the direct detection scattering cross section is organized as follows. First, we derive 
the non-relativistic interaction Hamiltonian for dark atoms and standard model nucleons from Eq. (j33p . 
Second, we use this to calculate the differential cross section for a dark atom to scatter from a spin singlet 
to a spin triplet state off of a standard model nucleon and append a form factor to account for recoil of the 
entire nucleus. Third, we rewrite the resulting rate in terms of the nuclear recoil Eji. Finally, we convolve 
the recoil rate with the dark matter velocity distribution. 

3.1 Non-relativistic Interaction Hamiltonian 

In order to calculate the scattering cross section for dark atoms off of standard model nuclei we derive the 
interaction Hamiltonian by taking the non-relativistic limit of the current-current interaction 

A =< J^D^^J^ >, (35) 
where D^^ is the Coulomb gauge propagator for X. The leading behavior of Eq. (f35l) is 



A 



+ 



ff5 e Ct» e 'I 

((7e • q){q-Pn) 

m„(g2 + Ml) 



2 me 



2 m„ m. 



• q){q-Pn 



+ 



o"p • q 

2m„ 



(<?p • q)q' 



2m„(g2+M2.) m„(g2 + M2.) 2m„(g2 + M; 



Xs ^r' 



(36) 



where m„ is the nucleon mass, pe is the initial momentum of the dark electron, that of the nucleon is pn 
and g is the momentum conjugate to the relative coordinate between the electron and nucleon. Xs^ndx^' 
are the initial and final spin states of the atom and can be written in the form: XAtom = Xp Xe- The 
dark matter spin operators are 5e,p = 1 <7e/2, c?p/2 (8) 1. £,r and S,r' are the initial and final spin states 
of the standard model nucleon. In the following analysis we consider proton masses of 0(100 GeV) and 
electron masses of O(lGeV), so we have ignored terms suppressed by rup. We have also dropped terms 
suppressed by Mz- Finally, we have omitted terms which depend on the spin of the standard model 
nucleon, as we expect these terms to contribute incoherently to the overall scattering cross section and 
hence be suppressed by the atomic number of the nucleus. 



3.2 Inelastic Dark Atom - Nucleus Scattering 

The cross section in the center of mass of the hydrogen-nucleus system is given by 



dn ^ dn^ 



0^ S = l) 



47r2 



(pii, p^, A^', H'\ H,nt Iph, P^, H) 



(37) 



Here Hint is the interaction Hamiltonian, which is obtained from Eq. (j36p by Fourier transforming q to 
position space. The prefactor contains una = ('^H + AfAf)/("iH + M]\f) which is the reduced mass of 
hydrogen and the nucleus, and the momenta k and k' which are the initial and final momenta conjugate 
to the relative coordinate between the atom and the nucleus. Our basis is the hydrogen atom momentum 
{Ph ' Ph}' nucleus momentum, {p^v i PatIj ^-^d the internal states of hydrogen and the nucleus, 
{A , A'} and {A^ , A^'}. The explicit evulation of the matrix element in Eq. (j37p is complicated by the fact 



9 



that there are four particles in the incoming and outgoing states. The hydrogen atom and the nucleus are 
both free particles while the electron and nucleon are bound the respective free particle motion. Since the 
scattering centers, the electron and nucleon, do not correspond to the coordinates of free particle motion, 
hydrogen's center of mass and the nucleus's center of mass, the matrix element in Eq. (j37p contains an 
atomic form-factor in addition to the usual nuclear form factor. We find, ignoring terms suppressed by 
me/rrin and mn/Mf^f, the cross section to be 




Fel{q 





2 




Xs 







Feliq' 



(38) 



The electron form factor Ff.i{q^) is found to be 



Fel{q^ 



1 + 



-2 



(39) 



where |0) is the ground state of the dark atom. The function Fuiq^) is the Helm nuclear form factor 
which accounts for the overlap between nucleon and nuclear states |4H I42j. We have averaged over initial 
nucleon spin and summed over final nucleon spin. Summing over the final atomic spin states and using 
Eq. ([39]) we have 



dn 



95 Ze e Cu 



I 



vr 



Fei Fh G{q^ 



\ 



Ml 



1 + 



Ml 



q 



(40) 



we have defined the function G{q'^) as follows 



G{q^) = 1 + ^ 

rUn 



F-i' - (1 + Ff 



q^+Ml 



(41) 



We rewrite the above in the following form 



dEn file ftff vl,i (1 + 2MmEr/MIY 



G\Er) 



eff 



2 (ff5 ec^)^' 

where we have defined the scale feff for compactness and for comparison to Ref. [IBj . 



(42) 
(43) 
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3.3 Modulated Nuclear Recoil Rate 



The amplitude of the modulated recoil rate at a detector with Nt target nuclei of mass is given by 



dR 



P 



-Nt 

2 niDM Jv, 



dVrel Vrel f{Vrel)-rEr- 



June 



January 



(44) 



where Vrei is the relativity velocity between the dark atom and the nucleus, p = 0.3 GeV/cm^ is the local 
dark matter density and we use fesc — 550 km/s for the dark matter escape velocity ~ see Ref. [431 . The 
lower bound on the velocity integration is the minimum relative velocity that can produce a given recoil 
energy 



Following Ref's. 



1 



(Mn + M-r 



j, the velocity distribution in Eq. (j44p is approximated by 



(45) 



The normalization is 



) J iri^f'e-^'^^'^, for.< 



fiv 



Ve 



N = erf 



for V > Vf .r.. 



2_ Vesc ^-{vesc/vof 



TT Vo 



(46) 



(47) 



with vq = 220 km/s. 



Figure [2] is an example of the modulated count rate at DAMA as defined in Eq. (j44p with the data points 
and reported uncertainties from DAMA and DAMA/LIBRA [5]. We have plotted the modulated spectrum 
for three choices of the set of parameters [m-p, rrie, a^, Mx, 55 and e) which satisfy the rather stringent 
list of constraints enumerated below. Note the linear dependence on Er and the presence of an atomic 
form factor in Eq. (j42p . Although the first term tends to push the peak toward larger values of Er, the 
atomic form factor turns off scattering when qao ~ 1. 



Mixing With The Standard Model: Perhaps the harshest constraints are on the mass of the axial 
U{1) and its kinetic mixing with the standard model, since the direct detection cross section is roughly 
proportional to e^/Mj[-. The one loop contribution of X to the anomolous magnetic moment of the muon 
is 



X 



a e 



dz ■ 



mi 



{l-zf 



+ m z 



(48) 



As discussed in Ref. j45l , regardless of how one treats the hadronic contribution to the theoretical prediction 
of a^, the X boson's one loop contribution must satisfy 



< 7.4 X 10" 



(49) 
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Modulated DAMA Spectrum 



R(cpd/kg/lceVee) 
0.025 




Figure 2: Examples of the modulated spectrum at DAMA defined in Eq. (j44p complete with the data 
points from the DAMA and DAMA/LIBRA experiments. The curves correspond to the following choice 
of parameters from left to right in order of their rise from zero (black, red and green in color order) 
rup = 200,100,70 GeV; nie ^ 2.1,1.7,1.7 GeV; feff ^ 103,92,67 GeV. The hyperfine splittings are ah 
about 111 keV. 



In order to be conservative we restrict ourselves to (see Figure 1 in Ref. [45j ) 

Mx > 100 MeV and < 10~^ 

Mixing between the massless gauge boson in the dark sector and the photon is not induced by loops in 
our theory, and we nominally set it zero. A constraint on this mixing, e', can be derived from bounds on 
its contribution to the anomolous magnetic moment of the electron. The constraint is e' < 0(10"^) [45j. 
Astrophysical constraints also exist, but are much less restrictive for the range of electron masses we are 
considering [46j . 

Sufficient Recombination: The residual ionized dark matter will scatter elastically as it does not cost 
any energy to flip a free spin. Efficient recombination and a hyperfine splitting consistent with DAMA 
imply that the typical electron mass is 0(1 GeV) and therefore too small to induce observable nuclear 
recoils. The strongest constraints on direct detection of the free dark protons come from CDMS [6t I47j . 
With a net exposure of 174.7 kg-d, the CDMS experiment allows 5.3 signal events at 90% confidence level. 
To be consistent with the bounds from direct detectioiH we demand that < 10"'^. 

Energy Level Corrections Due to the Axial U(l): The proton-electron interaction due to the 
broken axial U(l) is a perturbation to the hydrogen Hamiltonian and gives a correction to the hyperfine 
level splitting. The correction is given by 

^The discussion in this section actually puts a bound on the local ionized fraction. We assume for simplicity - here and 
throughout the paper - that the distribution of ionized dark matter matches that of atomic dark matter. However, due to the 
presence of a long-ranged force, it may be that the ionized distribution is very different from the typical dark matter halo. A 
full N-body simulation of a multiple species halo is beyond the scope of the present work, so for now we ignore this interesting 
possibility. 



12 



Mp = 100 GeV Mp = 200 GeV 




Figure 3: The allowed parameter space for Mx and for two values of the dark proton mass. The other 
atomic parameters - an and rrie - have been chosen so that E^f 111 keV and < 10~^. On each plot, 
the displayed values of are 10~^, 10~^ and 10~^ from left to right (green, red and blue in color order). To 
be consistent with constraints discussed in Ref. 05], Mx > 0.5 GeV for = 10~^. The allowed points have 
an average rate, in the 2 to 6 keVee bins at DAMA, between 0.99 x 10^^ and 1.63 x 10^^ cpd/kg/keVee. 
The excluded regions correspond to choices of parameters which do not satisfy Eq. ([53]); the solid black 
line is for A = 1 TeV and the dashed line is for A = 10 TeV. 



6Ehf ~ (0| ^ 10) = /^(l + aoMx)-^. (50) 

47r r 47r oq 



Requiring 6Ehf <C -E/i/ gives 



2 32 7r 3 rrif. Mx 



2 



gi^—-ai^[l + ^ . (51) 

3 rup \ ODfJ-H, 



The parameter sets shown in Figure [2] satisfy this constraint. 

Breaking the Axial U(l): Masses for the dark electron, dark proton and the axial gauge boson all 
violate the axial U(l) symmetry. Perhaps the simplest way to give mass to these particles is giving a vev 
to a charge +2 scalar (j)+2, as in [IB] 



LI 2 / 2 2 \ ^ 

Cb\D^cI)+2\ -X{\cI)+2\ -v+2) +yp{PR4'+2PL+PL(l)*+2PR) +ye{eB^(l>+2eL + eL(l)*+2(^R)- (52) 
When the scalar is at its vev, the mass spectrum is 

Mx = 95 V+2 
rrip = Vp v+2 

rUe = yeV+2- (53) 

The DAMA signal requires nip > mx and 55 > O(10~^) thus one might worry about the perturbativity 
of Up. The yukawa coupling runs according to the following one loop renormalization group evolution [48] 
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which blows up at the scale A. If we take A = 1 TeV or 10 TeV, our parameter space is constrained as shown 
in Figure [3l In principle, the proton could be a composite object and the axial-symmetry breaking could 
occur at strong coupling (as in QCD) and not via a weakly coupled scalar. The proton could also carry a 
charge under another gauge interaction that is relatively strong, but breaks at a TeV, thus tempering the 
UV behavior of yp. We leave explicit models of UV completions to future work. 

Figure [3] displays the allowed parameter space for a few choices of Mx, 95 ande with < 10^*^ level. 

4 Discussion 

Dark matter succinctly explains a number of astrophysical and cosmological observations that are other- 
wise puzzling. Standard WIMP dark matter can accommodate the gross features of these observations and 
naturally exists in models that attempt to explain the origin of the weak-scale. However, the typical WIMP 
seems unable to explain observed small-scale structure and tensions between direct detection experiments. 
These considerations point to the possibility of a non-minimal dark sector, which contains more similarities 
to the light sector than is typically thought. Atomic dark matter - with a non-negligible ionized fraction 
Xe and a new massless gauge boson - offers the possibility of significantly different phenomena in the dark 
sector than those of standard WIMPs. 

As discussed in Section O the residual ionized fraction can keep dark matter in equilibrium with dark 
radiation long enough to smooth halo structure on small scales. Furthermore, atomic dark matter may 
have hyperfine transitions of the right size to offer an inelastic explanation for the DAMA data. However, 
having X^ large enough to smooth out structure is inconsistent with the simultaneous postitive DAMA 
signal and null CDMS signal under the assumption that the charged halo has the same distribution as 
the atomic halo. Yet, if X^ is large enough, the distribution of the ionized fraction may be smoother and 
more spherical than standard halo models as suggested in Ref. |49l I50j . If the local ionized dark matter 
distribution is very different from the atomic dark matter distribution, then parameter space exists which 
can explain both small-scale structure and the DAMA signal. On the other hand, if our simple model of 
atomic dark matter is the right explanation for DAMA and the ionized components of the halo follow the 
distribution of the atomic dark matter, then other direct detection experiments should see dark protons in 
the near future. Simulations of stucture formation with charged and neutral components could shed light 
on these issues. 

The dynamics that lead to atomic dark matter also may have other phenomenological implications. For 
example, in parts of parameter space where the ionized fraction is large enough, H2 molecules may form 
through processes catalyzed by the residual ions, as in the SM [51] 



H + e <-> H +7 
H + H ^ H2 + e 



(55) 



and 



H + p ^ H+ + 7 
H+ + H ^ H2 + P. 



(56) 
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The existence of molecular states in the dark sector offers the possibility of cooling mechanisms which, in 
the SM, are thought to be very important for the formation of the first stars [52]. This raises the interesting 
question of whether and to what extent compact objects, e.g. dark stars, could form for weak-scale dark 
atoms. Moreover, if the dark photon mixes with the SM photon, it may result in dark atomic line emissions 
in cosmic gamma rays. 

We have presented a somewhat generic model of atomic dark matter. Explicit models which explain the 
asymmetry abundance and which serve as ultraviolet completions of the model could potentially relate 
astrophysical phenomena to physics to be probed by the Large Hadron Collider. The part of parameter 
space in which the measured DAMA signal is post-dieted requires the dark proton to be strongly coupled, 
or nearly so, at a TeV. If strongly coupled, one could imagine additional features of the dark sector - i.e., 
a composite atomic nucleus - which more strongly mimic our visible world. 

The authors would like to thank Julian Krolik, Kirill Melnikov, Colin Norman, and Alex Szalay for helpful 
discussions. This work is supported in part by the National Science Foundation under grant NSF-PHY- 
0401513, the Department of Energy's OJI program under grant DE-FG02-03ER4127, and the OWC. 
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